ALZHEIMER'S DISEASE: ANALYSIS OF A MATHEMATICAL MODEL 
INCORPORATING THE ROLE OF PRIONS 



MOHAMED HELAL, ERWAN HINGANT, LAURENT PUJO-MENJOUET, AND GLENN E WEBB 



Abstract. We introduce a mathematical model of the in vivo progression of Alzheimer's disease with 



focus on the role of prions in memory impairment. Our model consists of differential equations that de- 
scribe the dynamic formation of (3-amyloid plaques based on the concentrations of A(3 oligomers, PrP*- 
^ proteins, and the Ap-x-PrP*- complex, which are hypothesized to be responsible for synaptic toxicity. We 

?H prove the well-posedness of the model and provided stability results for its unique equilibrium, when the 

polymerization rate of |3-amyloid is constant and also when it is described by a power law. 
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1. Introduction 

1.1. What is the link between Alzheimer disease and prion proteins? Alzheimer's disease (AD) 
is acknowledged as one of the most widespread diseases of age-related dementia with « 35.6 million 
people infected worldwide (World Alzheimer Report 2010 |40|). By the 2050's, this same report has 
predicted three or four times more people living with AD. AD affects memory, cognizance, behavior, 
and eventually leads to death. Apart from the social dysfunction of patients, another notable societal 
consequence of AD is its economic cost (w $422 billion in 2009 |,40J). The human and social impact 
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of AD has driven extensive research to understand its causes and to develop effective therapies. Among 
recent findings are the resuhs that imply cellular prion protein (PrP*^ ) is connected to memory impair- 
ment ||7][8l [r7{[23]|27) . This connection is the focus of our modeling here, which we hope will contribute 
to understanding the relation of AD to prions. 

The pathogenesis of AD is related to a gradual build-up of |3-amyloid (A(3) plaques in the brain 
1 13 20). (3-amyloid plaques are formed from the AP peptides obtained from the amyloid protein 
precursor (APP) protein cleaved at a displaced position. There exist different forms of (3-amyloids , 
from soluble monomers to insoluble fibrillar aggregates ||5]|25]|26j|37]|^l. It has been revealed that 
the toxicity depends on the size of these structures and recent evidence suggest that oligomers (small 
aggregates) play a key role in memory impairment rather than (3-amyloid plaques (larger aggregates) 
formed in the brain |35|. More specifically, A(3 oligomers cause memory impairment via synaptic 
toxicity onto neurons. This phenomenon seems to be induce d by a membrane receptor, and there is 
evidence that this rogue agent is the PrP*-- protein | [l7l[23|[28][32]|33| We note that this protein, when 
misfolded in a pathological form called PrP^'^ , is responsible for Creutzfeldt-Jacob disease. Indeed, 
it is believed that there is a high affinity between PrP*^ and A (3 oligomers, at least theoretically |16|. 
Moreover, the prion protein has also been identified as an APP regulator, which confirms that both are 
highly related [28j,38J. This discovery offers a new therapeutic target to recover memory in AD patients, 
or at least slow memory depletion |[6][T4[. 



1.2. What is our objective? Our objective here is to introduce and study a new in vivo model of AD 
evolution mediated by PrP*-- proteins. To the best of our knowledge, no model such as the one proposed 
here, has yet been advanced. There exist a variety of models specifically designed for Alzheimer's 
disease and their treatment, such as in |[2[3]|4][T0|[TT][l5][T8l[T9]|24||31^ Nevertheless, the prion 
protein has never been taken into account in the way we formulate here, and our model could helpful in 
designing new experiments and treatments.. 

This paper is organized as follows. We present the model in section|2] and provide a well-posedness 
result in the particular case that |3-amyloids are formed at a constant rate. In section [3] we provide a 
theoretical study of our model in a more general context with a power law rate of polymerization, i.e. 
the polymerization or build-up rate depends on (3-amyloid plaque size. 



2. The model 

2.1. A model for beta-amyloid formation with prions. The model deals with four different species. 
First, the concentration of A(3 oligomers consisting of aggregates of a few AP peptides; second, the 
concentration of the PrP*^ protein; third, the concentration of the complex formed from one AP oligomer 
binding onto one PrP*-- protein. These quantities are soluble and their concentration will be described in 
terms of ordinary differential equations. Fourth, we have the insoluble P-amyloid plaques described by 
a density according to their size x. This approach is standard in modeling prion proliferation phenomena 
(see for instance ||4[[T8][30J). Note that the size x is an abstract variable that could be the volume of the 
aggregate. Here, however, we view aggregates as fibrils that lengthen in one dimension. The size 
variable x thus belongs to the interval (xq, +°°), where xq > stands for a critical size below which the 
plaques cannot form. To summarize we denote, for x E (xq, +°°) and f > 0, 

■ f(t;x) > : the density of (3-amyloid plaques of size x at time f, 

■ u{t) > 0: the concentration of soluble A|3 oligomers (unbounded oligomers) at time f, 

■ pit) >0: the concentration of soluble cellular prion proteins PrP*^ at time f, 

■ b{t) > 0: the concentration of A(3-x-PrP*^ complex (bounded oligomers) at time t. 
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Figure 1. Schematic diagram of the evolution processes of (3-amyloid plaques, 
A(3 oligomers (bounded and unbounded), and PrP*^ in the model. 



Note that |3-amyloid plaques are formed from the clustering of A|3 oligomers. The rate of agglom- 
eration depends on the concentration of soluble oligomers and the structure of the amyloid which is 
linked to its size. It occurs in a mass action between plaques and oligomers at a nonnegative rate given 
by p{x), where x is the size of the plaque. This is the reason why the intentionally misused word "size" 
considered here (and described above) accounts for the mass of A (3 oligomers that form the polymer 
We assume indeed, that the mass of one oligomer is given by a "sufficiently small" parameter e > 0. 
Thus, the number of oUgomers in a plaque of mass jc > is jc/e which justifies our assumption that the 
size of plaques is a continuum. Moreover, amyloids have a critical size xq = £n > 0, where n e N* is the 
number of oligomers in the critical plaque size. The amyloids are prone to be damaged at a nonnegative 
rate /j, possibly dependent on the size x of the plaques. All the parameters for A(3 oligomers, PrP*^ , and 
(3-amyloid plaques, such as production, binding and degradation rates, are nonnegative and described 
in table O 

Then, writing evolution equations for these four quantities, we obtain 

^^f{x,t)+u{t)^[p{x)f{x,t)]=-ij{x)fix,t) on(xo,+°°)x (0,+°o), (1) 

1 /■+" 

u = Xu—yiiU — lup + ab — nN{u) u p{x)f{x,t)dx on(0,+oa), (2) 

e Jxo 

p = 'k.p-'Ypp-'iup + ab on(0,+oo), (3) 



b^zup-{a + 5}b on(0,+oo). (4) 
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Parameter/Variable 


Definition 


Unit 


t 


Time 


days 


X 


size of [3-amyloid plaques 




Xq 


Critical size of |3-amyloid plaques 




n 


Number of oligomers in a plaque of size xq 




e 


Mass of one oligomer 






Source of A|3 oligomers 


days^' 




Degradation rate of A (3 oligomers 


days^' 




Source of PrP'^ 


days^' 


yp 


Degradation rate of PrP*^ 


days^' 


T 


Binding rate of A(3 oligomers onto PrP*- 


days^' 


a 


Unbinding rate of Ap-x-PrP'- 


days^' 


8 


Degradation rate of Ap-x-PrP'' 


days^' 


P{^) 


Conversion rate of oligomers into a plaque 


(SAF/sq)-'*-days-' 




Degradation rate of a plaque 


days^' 



Table 1 . Parameter description of the model. * SAF/sq means Scrapie- Associated 
Fibrils per square unit and is explained in detail by Rubenstein et al. |34j (we consider 
plaques as being fibrils here). 



The term accounts for the formation rate of a new p -amyloid plaque with size xq from the A|3 oligo- 
mers. In order to balance this term, we add the boundary condition 

u{t)p{xo)f{xQ,t)^N{u{t)), t>0. (5) 

The integral in the right-hand side of equation (|2| is the total polymerization with parameters 1 /e, since 
djc/e counts the number of oligomers into a unit of length dx. Finally, the problem is completed with 
nonnegative initial data, a function /"' > and u'", p"\ b'" > 0, such that at time f = 

/(•,?=0)=r on(^o,+-), (6) 

and 

u{t = 0) = u'", p{t=0)=p'" and b{t^O)=b"'. (7) 
The above system ([T]|5]l involves two formal balance laws: the first one for prion proteins 

d 



V- rp) = Xp-'Ypp-?)b, 



and the second for A (3 oligomers 



d ( 1 /■+" \ 1 /'+°° 

— Nj + mH — / xfdx =\,,—y,,u — 't)b / xufdx. 

dt\ eJxo ' J e Jxo 

The total concentrations of both evolve in time according to the production and degradation rates. In 
figure [T] we give a schematic representation of these processes. 

2.2. An associated ODE system. In this section we investigate constant polymerization and degrada- 
tion rates, i.e, rates independent of the size of the plaque involved in the process. This first approach 
is biologically less realistic, but technically more tractable, yet still quite challenging for an analytical 
study of the problem. In section [3] the polymerization rate p will be taken more realistically as a power 
of X. Here we assume that p(x) := p and ij{x) := /j are positive constants. Moreover, without loss of 
generality, we let e = 1, which only requires a rescaling of the units in the equations. Then, we assume a 
pre-equilibrium hypothesis for the formation of (3-amyloid plaques, as formulated in |29| for filaments. 



ALZHEIMER'S DISEASE: ANALYSIS OF A MATHEMATICAL MODEL 



5 



by setting N{u) = au". The formation rate is given by a > and the number of oUgomers necessary to 
form a new plaque is an integer, « > 1. With these assumptions we are able to close the system ([T]|4]i 
with respect to (|5]) into a system of four differential equations. Indeed, integrating ([TJ over (xq, we 
get formally an equation over the quantity of amyloids at time f > 

A{t) = I f{x,t)dx. 



This method has already been used on the prion model in p8) . Now the problem reads, for f > 0, 

A = au"-ijA, (8) 



u = Xii —JuU — Xup + ob — anu" — puA, 
p = Xp- jpp - Tup + ob, 
b = lup- {<3 + b)b. 



(9) 
(10) 
(11) 



The mass of (3-amyloid plaques is given by M{t) = J^"" xf{x,t)dx which satisfies an equation (formal 
integration of ([T])) that can be solved independently, since 

M ^ nau" + puA - ijM . (12) 

Notice that initial conditions for A and M are given by A'" = J^°° f'"{x)dx and M'" — J^°° xf" {x)dx, 
while the initial conditions for u, p and b are unchanged. 

The next subsection is devoted to the analysis of the system (8pTi. 



2.3. Well-posedness and stability of the ODE system. We prove in the following proposition the 
positivity, existence, and uniqueness of a global solution to the system (8pTi with classical techniques 
from the theory of ordinary differential equations(pT|). 

Proposition 1 (Well-posedness). Assume A,,,, Xp, y„, jp, z, a, 5, p and fi are positive, and letn> 1 be an 
integer. For any {A'",u'",p'",b'") £ there exists a unique nonnegative bounded solution {A,u,p,b) 



to the system (S]|77) defined for all time t > 0, i.e, the solution A, u, p and b belong to C^{R+) and 
remains in the stable subset 



S= I {A,u,p,b) eM+ : nA + u + p + 2b <nA 




(13) 



with X — Xu+Xp and m = min{/j, y„ , Jp , 5}. Furthermore, let M{t — 0) = M'" > 0, and then there exist 
a unique nonnegative solution M to \i2) , defined for all time f > 0. 



Proof. LetF 



1-^ . 



be given by 



F{A,u,p,b) 



Fi 

\F4 



■.= au"—pA \ 
■ Xi, — Y„M — xup + ab — anu" — puA 
- Xp — jpp - lup + ab 
■.= zup-{a + b)b ) 

F is obviously and locally Lipschitz continuous on M^. Moreover, if {A,u,p,b) G K^, f i > when 
A = 0, F2>0 when u ~0, Ft, >Q when p = 0, and > when b = 0. Thus, the system is quasi-positive 
and the solution remains in . Finally, we remark that 

^ {nA + u + p + 2b) < X-m{nA + u + p + 2b), 
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with A, = A-,, +Xp and m — min{/j,Y„,Yp,5} > 0, and Gronwall's lemma ensures that 



nA{t) + u{t) + p{t) +2b{t) <nA" 



-2b" 



This proves the global existence of a unique nonnegative bounded solution {A,u,p, b). The claim for 
the mass M is straightforward. □ 

We next consider the existence of a steady state A„, m„, b^ and the asymptotic behavior of 
solutions to ([8|TT]i. It is easy to compute the steady state by solving the problem 



- au'l = 

A-i, — JiiUoo — TMooPoo + oboa — GLnu'l, — pWooAoo — 

Xp — JpPoo — TMco/^oo + O/Joo = 

TMooPoo - (8 + a)fooo = 



(14) 
(15) 
(16) 
(17) 



From the structure of the second equation, we cannot give an explicit formula for this problem. To 
obtain Moo we have to solve an algebraic equation, which involves a polynomial of degree n. However, 
we can prove that the solution exists, and then Uoo is given implicitly. The next proposition establishes 
the local stability of the steady state.. 

Theorem 2 (Linear stability). Under hypothesis of proposition^ there exists a unique positive steady 
state Aca, Moo, Poo and boo to (S][77i with 



a 



1 Ap(T-T*) 



fj i:*Moo + Yp OT*Moo + Yp 

where 1* = t(1— o/(5 + o) and Moo is the unique positive root of Q, defined by 

Q{x) — ypXu +ax — P{x), for every x>0 
with a — 1* (A-j, — Ap) — YmYp '^"'^ 



.n+l 



P{x) ^ i*'YuX + ajpnx" + {ax*'n + pYp-)x' 
Moreover, this equilibrium is locally linearly asymptotically stable. 



-px —X 



,n+2 



Proof. First, equation ( [T4] l gives Aoo with respect to Moo. Then, combining ([T6| and ( [T7] l we get poo and boo 
as functions of Moo. Now replacing poo and boo in ([T5| we get Moo as the root of Q. It is straightforward that 
Q has a unique positive root. Indeed, it is the intersection between a line and a monotonic polynomial 
on the half plane. Now, we linearize the system in Aoo, Moo, poo and boo. Let X ~ {A,u,p,bY and the 
linearized system reads 

d 



-X=DX, 



where 



D = 



dt 

—/J anu'^^ 
-pwoo Ym - '^P«- - OLn^W^^ - pAo. 
-zpoo 

T/?oo 





-tmo, 



(Yp+tMo. 



IP 

xu, 




a 
a 

-(a + 8) J 



The characteristic polynomial is of the form 

P(A) = A'* + fliA^+a2A^ + fl3A + fl4, 
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with the fl,- > 0, / = 1 ... 4 given in the appendix. Moreover it satisfies 

2 2 

Then, according to the Routh-Hurwitz criterion (see |2] Th. 4.4, page 150]), all the roots of the charac- 
terisic polynomial P are negative or have negative real part, thus the equilibrium is locally asymptotically 
stable. □ 



To go further, we give a conditional global stability result when no nucleation is considered, i.e., 

a = 0. 

Proposition 3 (Global stability). Assume that a = 0. Under the condition 



the unique equilibrium is given by 



. „ \, , 1 A-p(t-T*) 

Aoo=0, noo = , boo = Moo, 

i*u„+yj, az*u^ + jp 

where Uoo is the unique positive root ofQ[x) = YpA,„ + ax — T*A-„x^, with a = X* (A-„ — A,p) — YuYp. Further, 
this equilibrium is globally asymptotically stable in the stable subset S defined in ( |13| l. 

Proof. The proof is given by a Lyapunov function <I> stated in the appendix. It is positive when the 
condition above is fulfilled and its derivative along the solution to the system ([8pT[) is negative definite. 
Thus, from the LaSalle's invariance principle, we get that under these hypotheses the equilibrium of 
(SfTTi is globally asymptotically stable. □ 



3. A POWER LAW POLYMERIZATION RATE 



The assumption that the polymerization rate p and the degradation rate fj are constant is not always 
biologically realistic, as recognized in |3j[15J. Consequently, we study here the more realistic case 
p(ji:) ^ x^, and in the following we restrict our analysis to e (0, 1). We will see that we are able to 
obtain a result of existence and uniqueness of solutions for this more general case. 



3. 1. Hypotheses and main result. We are interested in nonnegative solutions to the system ([T]|4]i with 
the boundary condition (|5]l, completed by initial data (|6| and (|7|, but with the new assumption p(x) ^ x^. 
Moreover, we require that our solution preserves the total mass of (3-amyloid in order to be biologically 
relevant. Hence, the solution / will be sought in the natural space L' (xq, +°°;xdx), since xdx measures 
the mass at any time. Our hypotheses for the system (fTpl are 



(HI) 



(H2) 
(H3) 
(H4) 



e (xo, +oo;xdx), f" > 0, a.e. x > xq. 



p > 0, p e w2-"([xo,-)), fJ>0, iJe wi-~([xo,-)). 



> , e W, ;"(M+), A^(0) = 0. 
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We note that (H2) implies the existence of a constant C > such that p{x) < Cx, with for example, 

C = 2||p'||l«' + p(xo) /xq. For any x > xq, we have 



PW < ||p'||l"(^ + xo)+p(xo) < 2||p'||z,' 



We remark that this kind of regularity of the rate p covers the case that p{x) ^ x^ with 9 G (0; 1). Also, 
(H3) implies the existence of a constant Km > such that N{w) < Kmw, for any w e [0,M]. Further, 
The nonnegativity of the parameters of table |2. 1 1 (hypothesis (H4)) is a natural assumption with regard 
to their biological meaning. 

We introduce the definition of a solution to system H^M- 



Definition 1. Consider a function f" satisfying (HI ) and let u'", p'", b'" be three nonnegative real data. 
Assume that p, /j, and all the parameters of table 2.1 verify assumptions (H2) - (H4), and let T > 0. 
Then a quadruplet {f,u,p,b) of nonnegative functions is said to be a solution on the interval {0,T) to 
the system with the boundary condition (|5]l and the initial data ^ and (|7|, // it satisfies, for any 
cpe CdO'T'] X [jco,+°°)) andte {0,T) 



r 



f{x,t)(p{x,t)dx - 



f"{x)(^{x,0)dx^ 
f{x,s) 



Jxo 



N{u{s))^){xo,s)ds 

^^(^{X,S) + U{s)p{x)^(^{x,s) -^(x)(p(x,5) 



dxds, 



and 



u{t) = m'" + 



p{t) 

b{t) 



-P 

: b'" 



— JiiU — lup + ob — xqN{u) — u 
[Xp — Jpp ~'lup + Gb] ds, 
zup— {G + 5)b] ds, 



p{x)f{x,s)dx 



ds, 



with the regularity f £ LT (O, T\L} (xo, +°°\xdx)^ and u,p,b £ ^"(0, T). 

Tlieorem 4 (Well-posedness). Let f" be a nonnegative function satisfying (HI), let u"', p" 
be nonnegative real numbers, and assume hypothesis (H2) to (H4). Let T > 0. There exists a unique 
nonnegative solution (/, u,p,b) to ([T]|4| with (|5]) and initial conditions given by (|6]l and (|7]l, in the sense 
of definition^ such that f £ CP {jO,T],L^ (xq, +°°;x''dx)^ for every r £ [0, 1], and u,p,b £ Cl{0, T). 

The proof of the theorem|4]is decomposed into two parts. First, we study the initial boundary value 
problem 

|;/(x,f)+M(f)|-[p(x)/(x,f)] ^ -lj{x)f{x,t) on {Xo,+oo) X (0,+°o), 

u{t)p{xo)f{xo,t) N{u{t)), on (0,+°o), 

/(•,f-0)-/'", on(xo,+oo). 
We prove in the subsection |3 .2 1 the following proposition: 



(18) 
(19) 
(20) 



Proposition 5. Let u £ C^(M+), let f" satisfy (HI), and assume hypothesis (H2) to (H3). For any 
T > 0, there exists a unique nonnegative solution f to ( |7(S|20| > in the sense of distributions, such that 
/ e ([0, r],Ll (jco, +oo;x''dx)) for every r £ [0, 1]. 
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The proof is in the spirit of the proof in f9) for the Lifshitz-Slyozov equation. It consists of a proof 
based on the concept of a mild solution in the sense of distributions, with the additional requirement of 
continuity from time into L' {xdx) space. 

The second step of the proof of theorem [4] is performed in subsection 3.3 Precisely, once we have 
the existence of a unique density /, when u is given, we are able to construct the operator 

S : C»([0,r])3 ^ C^i[0,T])^ 

(21) 

(u,p,b) {Su,Sp,Sh) ^S{u,p,b), 



Su 

Sn 



= P 



Sh = b" 



K, — juU — xup + ab — xqN{u) 
\kp — Jpp — 'lup + ob] ds, 
[lup — {o + 5)b] ds, 



p{x)f{x,s)dx 



ds, 



where / is the unique solution associated to u given by proposition |5] Then, theorem|4]is finally proven 
in subsection 3.3 applying the Banach fixed point theorem to the operator S. 

3.2. Existence of a solution to the autonomous problem. In the following we let u g C^(M+) and 

we use the notations a{x,t) ~ u{t)p{x) and c{x,t) — ~u{t)p'{x) for every {x,t) G [xo,+°o) x R^. From 
(H2) and noting that p (x) < Cx, we have for any f > 



a{t,x) <Ax, forjic>xo, 
\a{t,x) — a{t,y)\ <A\x — y\, forx,y>xo, 
\c(t,x)\<B, 



.xo,- 



(22) 
(23) 
(24) 

In order to establish the mild 



where A = max(C|lM|lL^,||M|lL~||p'||L-) and B = |lM|lL»|lp'||i^ 
formulation of the problem, we define the characteristic which reaches x > xq at time f > 0, that is, the 
solution to 

— X(s;x,t) — a(t,X{s;x,t)), 

(25) 

X{t;x,t) =x. 

From property ( |23| ), their exists a unique characteristic that reaches (x,f ).We note that it makes sense as 
long as X{s;x,t) > xq. Thus, we define the starting time of the characteristic as 

so{x,t) infji e [0,?] : X{s;x,t) > xq} . 

The characteristic will be defined for any time s > sq and takes its origin from the initial or the boundary 
condition, respectively, if io = or sq > 0. We recall the classical properties of these characteristics 

X{s;X{a;x,t),a) ^X{s;x,t) 
d 



J{s;x,t) :— —X(s;x,t) = exp fj^ c{a,X{o;x,t))do 



-X{s;x,t) 



-a{t,x)J{s;x,t). 



Also, remarking that so{X{t;xo,Q),t) — 0, then by monotonicity and continuity ofX for any f > 0, we get 

X e (xo,X{t;xQ,0)) 4=^ S(){x,t) S (0,f), and for any x £ {xo,X{t;xo,Q)) we have X{so{x,t);x,t) = xq. 
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It follows that for every x e {xQ,X{t;xQ,0)) 

I{x,t) := -^soix,t) =J{so{x,t);x,t)/a{so{x,t),XQ). 
ax 

Considering the derivative of f{s,X{s;x,t)) in s, and integrating over {so,t) we obtain the mild formu- 
lation of the problem. The mild solution is defined for a.e. {x,t) E {xq, +°°) x by 



/"'(X(0;x,f))y(0;x,f)exp - / ij{X{o;x,t))da ] x>X{t;xo,0), 

f{x,t)^{ A., \ (26) 

N{u{sQ{x,t)))I{x,t)exp ( - / fj{X{G;x,t))da I x e {xQ,X{t;xo,0)). 

We infer from the formulation (|26| that for a.e {x,t) e [xo,+oo) x M+, / is nonnegative, since J and / 
are nonnegative, and /'" satisfies (HI). We recall some useful properties that are derived in |9, Lemma 
1]- 

Lemma 6. Let u € be a given data and assume that (H2) holds. Then for any x>xq and f > 0, 

as long as the characteristic curve s ^ X{s;x,t) defined in \25) exists, i.e., s > so(xj), we have 



for SI <S2, X{si;x,t) < X{s2;x,t) < X{si;x,t)e^^''^--''^'> 
if x„ — >■ +00, then for all t > s >0, X{s;xj) — > +°° 
fors>t, X{s;x,t) <xe'^^'''K 

Proof. We refer to |[9] Lemma 1], where the result follows from the fact that for any x > xq, t > and 

so(x,t) < 51 < 52, we have 



l'S2 rS2 

XQ ^ X(s2',x,t) — X{si;x,t) + / a{s,X{s;x,t))ds <X{si;x,t) +A / X{s;x,t)ds, 

Jsi Jsi 

where A is given by (|22]l. □ 



In the sequel we will repeatedly refer to the changes of variables 

y — X{Q;x,t) ower X G {X{t,XQ,Q),+°o), with Jacobian 7(0;x,f), 
s = so{x,t) over X E {xo,X{t;xo,Q)), with Jacobian —I{x,t). 
The first is a - diffeomorphism from {X (f , xq , 0) , +°°) into (xq , +°°), and the second from (xq , X (f ; xq , 0) ) 
into (0,f). Integrating / defined by (|26| over {0,R) with R > X{t;xo,Q), using the change of variables 
above, using lemma [6j and taking the limit R we get 



x\f{t,x)\dx< r X{t;y,0)\f'"iy)\dy+ f X{t;s,xo)\N{u{s))\ds 

JXQ Jo 

< e^' (^l%\f"iy)\dy + l\o\Niuis))\ds^ , 



(27) 



where we have split the integral into two parts and uses both the previous changes of variables. Thus,for 
any T > 0, f gL" (p,T;L^ (xq, +°o;xiix)), and therefore inL°° (O, T;L^ (xq, +°°;x''dx)^, for any re [0, 1]. 
In the next lemma we claim that / defined by (|26]l is a weak solution. 



Lemma 7. Let f be the mild solution defined by \26\ . Then for any t > 

f{x,t)(\){x,t)dx= [ /"'(x)(p(x,0)iix+ / N{u{s))(p{xQ,s)ds 
Jxn ' Jo 

f{x,s 



JxQ 



^^(p(x,5)M(5)p(x)^(p(x,5)-jU(x)(p(x,i) 



dxds, 



forall(f,eCr{[0,T]x[xo,+^)). 
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Proof. Since / belongs to / e L°° (O, T\0 {xQ,-\-°°;xdx)), it is possible to multiply the mild solution / 
against a test function (p e C^{[Q, T] x [xq, +00)) and integrate over (xq, to obtain 

°/(x,f)cp(x,f)t/x= f"{y)^,iX{t-,y,0))e-^oM{x(a-y,o))dc^y 

JXQ 

- /'A?(M(i))(p(X(f;xo,5),f)e-^'"(^(°'"«-^»'"'£/i, (28) 
Jo 



by the same change of variable made above for ( |27| i. Furthermore, we have 

rt rX(s-^xofl) 

/ / f{x,s)[dt(p{x,s) +a{s,x)dx(?{x,s) — ij{x)(p(x,s)]dxds 

Jo Jxn 



'-to 

rt /.+00 ^ 
.0 



(x)(p(X(r;x, 0) , f )e- /of (^(<';>'.o))''<^ _ Z^" /'« (jc)(p(x, 0)t/x, (29) 



'.to Jxo 

stiU using the change of variable mentioned above and 
rt 

f{x,s) P((p(x,5) +a{s,x)dx(f{x,s) — /j(x)(p(x,i)] c/xc/i 



/O JX(.s;xo.O) 



dzds 



f fN{u{z))^ ((p(X(5;xo,z),5)e-^'''(^("^^«'^))^-^ 

iO "5 V J 

f Niuis))(i?{X{t;xQ,s),t)e-^'^''^^^''-'<''''>'>'"'dzds- f Niu{s))(!;>{xo,s)ds. (30) 
Jo 



Finally, combining ( [28| ), ( [29| l and ( (30| l we obtain that / is a weak solution. □ 
The aim of the following lemma is to prove that the moments of / less than 1 are continuous in time. 



Lemma 8. Let hypothesis (HI) to (H3) hold. Let f be the mild solution given by \26\ . Then for any 
T>0, 

/eC°([0,r],L'(xo,+oo;x'c/x)), for every re [0,1]. 

Proof Let T > and r e [0, 1], since / e L^^^. (M+,L' (xq, +°°;x''dx)^, we have for any t > and 5f > 
such that t + St <T 



[ x' |/(x, f + 5f ) - /(x, f ) I c/x = /i + /2 + /3 , 

J xn 



where 

/l = 



/ x''|/(x,f + 5f)-/(x,f)|t/x, 

/■J>!'(r+8f;,To,0) 

/2=/ x''\f{x,t'Vht)~ f(x,t)\dx, 

Jxinxofl) 



x''\f{x,t + ht)- f{x,t)\dx. 

Jx(t+5t;xofl) 

Our goal is to prove that each term goes to zero when 5? goes to zero. We first bound 73, which results 
from the initial condition, since for x > X(f + 5f ;xo,0) > X(f;xo,0), it follows that 

.f (X(0;x, f + 5f ))7(0;x, r + 5f)e- 

(X (0;x, f )) |7(0;x, f )e- IoK{x(o,xj))da 



h = 



12 M. HELAL, E. HINGANT, L. PUJO-MENJOUET & G. F. WEBB 

Let fl" G with compact support supp{f^) C (0,/?e) and converge in L^{[xQ,+o°),xdx) to /'". We 
write It, as follows 

h^l\+ll+ll (31) 

where 



l\ = / x''|/'"(X(0;x,f + 5f))-/^"(X(0;x,f + 5f))| 

ix(f+8r;.vo,0) 

ll = ^^/"(^(O;^,? + 5f)y(0;;c,f + 5f) 

ix(r+8?;xn.O) 



/X(r+8?;xo.O) 

X g-/o'+^'A-(X(a;x,r+&))rfa 

-/^"(X(0;x,f))y(0;x,f)e-i'o^(^('''^''»'"'|t/;c, 
ll = r x'-|/r(X(0;x,f))-.r(X(0;x,f))| 

Jx(f+8f;xo.O) 



Dropping the exponential term, which is bounded by one, and changing of variables y — X(0;x,f + 5f ) 
in and y — X{0\x^t) in /|, we get 



/' +il < 2e^^ r'°f\f{y)-ft{y)\dy = c^(r,8), (32) 

with the help of lemma |6] Next we bound by 

ll < r x''\fi"{X{0;x,t + dt))-fi"{X{0;x,t))\J{0;x,t + ?,t)dx 

Jx{t+&t;xo.O) 

x''fi"{X{0-x,t))\J{0;x,t + 5t)~J{0-x,t)\dx 

X{t+bt;xofl) 

x'fl"{X{Q;x,t))J{Q;x,t) 

X(f+Sf;xo,0) 

X |g-/o+''V(Jf(o;.x-,r+8r))rfo _ g-/^M^(o;x,OVa|^^^ 

and we denote the integrals by j\ to y|, respectively. We remark that 7(0, x, t) < e^^ by (|24]) and so 

-/g <e^^||/£"llz.- / x''\XiO;x,t + 5t) - X{0;x,t)\dx 

Jx{t+&t;xo.O) 



<dte''^\\fi"\\L- r sup 

ix(r+5f;.vo,0) sG[r,f+8f] 



|x(0;x,.) 



<htAe^^'^\\f^'\\L- r\'+^dx, (33) 



where Cg depends on T, A and 7?g /.e., the compact support of f". Then 

^3 < e^^ll/i'lk- xHe^''^^''^' - 1 \dx 

Jx(t+ht;XQfi) 
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with 



rt+bt rt 

|G(f,5f,x)|= / c(o,X(o;x,r + 5f))iio- / c{<3,X{<3;x,t))d(s 
Jo Jo 

rt+ht 

< / p'{X{a;x,t + dt))-p'{X{a;x,t)) u{a)da 
Jo 

rt+bt 

+ c{a,X{a;x,t)) 



Thus, with d22l and d24^ 



|G(f,5f,x)| < /r||M||i- / X{a;x,t + 5t)-X{a;x,t 
Jo 

< dtK\\u\\L-^ / sup 

■^0 s£[t,t+8t] 

<?>t{K\\u\\L-'ATe'^^x + B) , 



dt 



X{<3;x,s) 



do + 5tB 
da + 5tB 



where K is the Lipschitz constant of p'. Since x < R^, let CdT, e) = K\\u\\i«.ATe^^R^+B, and if |x| < y, 
then 



Thus, we get 



Since /j is nonnegative, < 



X{t+5t;xo.O) 



5i;'"p{X[0-^,t+U))d0-!'Qf,{X(G-^l))dG) 



dx. 



Exactly as above, 

rt+ht rt 

/ ij{X{o;x,t + 5t))da- ij{X{o;x,t))do <8tMATe'^^x + 5t\\ij\\L-, 
Jo Jo 

with M = Lipschitz constant of /j. Denoting by Cm{T,e) = MATe^^R^ + we get 

V / JxQ 

From ( (32] l, ( (33] l, ( (34| l and ( |35| l we can conclude that for any e > 0, 

h{ht)<C\{T,^)+Cl{T,?>t,z), 
with lim£^oC^(r,e) = and lim5,^oC|(r,5f,e) = 0. 



(34) 



(35) 



(36) 
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Next, concerning Ii , / can be written from the boundary condition. Let e such that — > u 
uniformly on [0, T]. Then we write I\ as follows: 

h< I x''\N{u{sQ{x,t + ht))~N{u^{sQ{x,t + ht))\I{x,t + ht)dx 

+ x'^ N{u^{so{x,t + 5t))I{x,t + 5t)e J-'o('-'+50''^ " 

Jxn 



/\N{u{s(j{x,t))-'N{u^{s{){x,t))\I{x,t)dx. 



dx 



From (H3) we obtain, similarly to It,, that there exist two constants C}(r,e) and C\{T,bt,€) such that 

/i(5r)<c;(r,e) + C?(r,&,e), (37) 

with limE^oC}(r,e) = and lim8,^oC^(r,5f,e) = 0. 
Finally, for h, we use the two formulas of /, 



h = 



/ x' N{u{soix,t + 5t)))Iix,t + 5t)e ^■'o("+5f)''^ ^ •' " 

Jx{t;xo,Q) 



dx 



Using the Lipschitz constant of denoted by K^, from the definition of / and with the help of lemma 
|6] we get 



h < x'oe'^'^+'^^'^KN |X(f + 5f;xo,0) -X(f;xo,0)| 



Using the regularization of /'", there exist two constants C2{T,e) and C2{T,dt,e) such that for any 



rX{t+5t;xo.O) 

+ x5e'''*^ / \f'"{X{Q;x,t))J{Q;x,t)\dx. 
Jx(f^ofi) 



e>0, 

/2(5f)<Ci(r,e) + cf(r,&,e), (38) 
withlimE^oC2(r,8) =Oand \\m^,^QCl{T,ht,E) =0. 

In conclusion, combining (l36|, ([37]i and (|38]l, we get for any £ > and 8f > 0, 



r x''\f{x,t + ht)^f{x,t)\dx<C\T,E)+&{T,ht,E), 

JxQ 

where C' (r,e) andC^(r, 5f,e) are two constants such that lim£^oC {T,^} = and lim5,^oC^(r, 5f,e) = 
0. Noticing that the proof remains the same when 5f is negative, taking the limsup in bt we get 

0<limsup f x''\f{x,t + 5t)~'f{x,t)\dx<C\T,E), for any e > 0. 

The proof is completed by taking the limit as £ goes to zero, which yields to the require regularity, 

/eC°([0,r],Li([xo,+°°),x'"t/r)forallre [0,1]. □ 

We finish this section with a useful estimate for the uniqueness investigation. 

Proposition 9. Let T > and u\, U2 G Cj^{Q,T). Let f\ and f2 be two mild solutions to ([T8|)-(|20]i, 



associated, respectively to ui and U2, with initial data /f", /" given by formula ( |26| l. Then, for any 

fG(o,r) 
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r'^x i/i (x, f ) - /2 (x, t)\dx< r°°x \fi" (x) - f'"{x) I dx 



Jxo 



lj{x)x (jc, i) — (x, i) I dxds + A 1 



7x0 



:\fi{x,s)-f2{x,s)\dxds 



+ (/^1,2+C||/2(vOIIl1(..</.))|«iW-«2WM5, 



where Ai is given by ^22\ for ui andKi2 is the Lipschitz constant of N on [0,R] withR = ma.x{\\ui\\[j^(^Qj^, 
II"2||l°°(0,7'))- Finally C > denotes a constant such that p{x) < Cx. 

Proof. This estimation is obtained from a classical argument of approximation. Let h — fi—f2 and 



h{x,t)(p{x,t)dx : 



h'"{x)(^{x,Q)dx+ / {N{ui{s)) - N{u2{s)))(?{xo,s)ds 



t /■+=< 



Jxo 



h{x,s) 



^^(p{x,s) + ai{s,x) ^(l^{x,s) ~ fi{x)(l?{x,s) 



dxds 



T 

Jxo 



(fli {s,x) — a2{s,x)) f2{x,s) —^){x,s)dxds. 

ax 



Let /!£ be a regularization of and 5g a regularization of the function. Take (p(x,i) — S^{h£{s,x))g{x) 
with g e C"([xo, +°°)). Then, letting 5 — > and then e — > 0, we get 



r \h{x,t)\g{x)dx = r \h"\x)\g{x)dx 

Jxo Jxo 



\N{uiis)) - N{u2{s))) Sign{ho{xo))g{xo)ds 



7x0 



\h{x,s)\ 



ai{s,x)—g{x)-/j{x)g{x) 
ax 



dxds 



(ai {s,x) ~ a2{s,x)) f2{x,s)Sign{h{s,x)) ^g{x)dxds. 
Jxo ox 

Finally, we approximate the identity function with a regularized function rj^ e C"([-^o, +°°)) such that 
r|^ (x) = X over (0,7?), and then taking the limit R — > +oo ends the proof. 

□ 



We get straightforward from proposition |7] that / defined by (|26| is a weak solution and the only 
one from propositio n [9| Indeed, getting ui — U2 and /f = f^ in proposition |9] leads to the uniqueness. 
Finally, proposition I8J provides the continuity in time of the moments with order less or equal to one. 
This concludes the proof of proposition |5] 

3.3. Proof of the well-posedness. In this section we prove theorem |4j We first study the operator S in 

Lemma 10. Consider hypothesis (H2) to (H4). Let u'", p'" and b'" be nonnegative initial data, and let 
f" satisfy (HI). Let M > be large enough such that u'",p'",b'" < M/2 and define 

Xm = {{u,p,b) e C°{[OJ])^ : < u,p,b<M} 

where C^{[0, 7"])^ is equipped with the uniform norm. Then, there exists T > (small enough) such that 
S : Xm ^ Xm is a contraction. 



16 



M. HELAL, E. HINGANT, L. PUJO-MENJOUET & G. F WEBB 



Proof. Let M be sufficiently large such that max{u"^,p"^,b"') < M/2, and let T > be small enough 
such that 

(Y„ + tM + o+xqCi (M) + C2(M, T))MT < M/2, 

{'Yp + zM)MT <M/2, 

{a + 5)MT <M/2, 

{K+oM)T <M/2, 

{Xp+aM)T <M/2, 

xM^T < M/2, 

where Ci (M) is the Lipschitz constant of on (0,M) and 

C2(M, r) = Ce''^^ (||.riLi(,^,) +Ci(M)Mr) , (39) 

where C is the constant such that p{x) < Cx, see p7| ). This assumption ensures that for any {u,p,b) £ 
Xm, then S{u,p,b) e Xm, i-e, the solution is bounded by M and is nonnegative. It remains to prove that 
5 is a contraction. Let {ui,pi,bi) and {u2,P2,b2) belong to Xm- Then 



I -^M I '^!( 



<'YuT\\ui -M2||<>o + Tr||Mi;?i -M2;'2lU+or||/7i - 

+ XorCl(M)||Ml-M2||oc 



-b2\ 



-T sup 



ul p{x)fi{x,s)dx~ u2 / p{x) f2{x,s)dx 

Jxn Jxn 



Then, 



sup 

ie[o.T] 



\\uiP\-U2P2\\oo <M\\ui-U2\U+M\\pi- p2\ 



(40) 



(41) 



m1 



/ p{x)fi{x,s)dx — u2 I p{x)f2{x,s)dx 

J XO J XQ 



<C2{MJ)\\ui-U2\\^ + CM sup 

(£[0,7-] 



x\fi{x,t) - f2{x,t)\dx 



and from Proposition |9j 



sup 



'XQ 



x\f\{x,t) - f2{x,t)\dx 



<T{Cx{M)+CC2{M,T))\\ui-U2\ 



(42) 



(43) 



We get similar bounds for \Sp^ — Sp^ |oo and \St,^ — |oc. We infer that there exists a constant C(M, T) 
depending only on M and T such that 



\{S„i,Spi,Sbi)- (5„2,5p,,5fo2)||oc < C{M,T)T\\{ui,pi,bi) - (m2,/'2,^2)I1oc, 



(44) 



with C{M, T)T ^0, when T goes to 0. Hence, if T is small enough such that C(M, T)T <l , then S 
a contraction. 



From Lemma 10 we have a local nonnegative solution on [0,T], which is unique with the solution 
{u,p,b) bounded by the constant M. The solution satisfies / e C''(0, T,L^ {xdx)) and u,p,b G C^{0, T). 
Futhermore from (H3), is continuous and from (H2), p{x) < Cx where C is a positive constant. Thus 
pf e C^{0, T;L} (dx)). We conclude that u,p and b defined in definition [ijhave continuous derivatives. 

Now we remark that the solutions satisfies on [0, T] 
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^ (m + /9 + 2b) = X„ + Xp -JuU- Jpp - 52/7 - nN{u) 



1 /'+°° 

u p{x)f{x,t)dx <'k — m{u+ p + 2b), 



e 



with m = min(Y„,Yp,5) and A- = A„ + Xp. Using Gronwall's lemma, the solutions remain bounded, at 
any time by 

u + p + 2b< u'" + p'" + 2b'" + - . (45) 

m 

From this global bound on u, p and b, we can construct the solution on any interval of time by repetition 
of the local argument. The proof of the theorem is complete. 



4. Summary 

The connection of prions and AD is not fully understood, but recent research suggests that soluble 
Ap oligomers are possible inducers of AD neuropathology. The key element of this hypothesis is that 
p— amyloid plaques increase their size over disease progression by the clustering of Ap oligomers, 
which are bound to PrP'-' proteins. Ap oligomers exist both as bounded and unbounded to PrP*-' proteins, 
and the agglomerartion rate in the formation of p— amyloid plaques depends on the concentrations of the 
bound and unbound Ap oligomers, the concentration of soluble PrP*^ , and the size of the p— amyloid 
plaques. We have introduced a mathematical model of the evolution of AD based on these hypotheses, 
and presented a mathematical analysis of its fundamental properties. Specifically, we have analyzed 
in detail the existence and uniqueness properties of solutions, as well as the qualitative properties of 
solution behavior. In specific cases we have quantified the stabilization of the solutions to steady state, 
a well-known feature of AD progression. In future work we will explore applications of this model to 
specific AD laboratory and clinical data. 

Appendix A. Characteristic polynomials of the linearized ODE system 

Here we give the coefficient a,-, / = 1 , . . . ,4 for the characteristic polynomial of the linearized system 
in proposition |2] 

V T*Moo + Yp A* J 



fl2= (A' + Y» + «n «~ (Yp + i:M=<. + o + 5)+Ypa+(Yp + TMoc)5 

^fjlju+T:--^, — +an^u'L-^ +P-u'l] +panul, + T:{yp + d)--^^ — , 
fl3= ( /J + Y» + an^Mr' +p-M» ) (YpO + (Yp + i:moo)5) + (Yp5 + (Yp + 5)a()t- 



;U ( Y, + o.n^u"J ' + p j + panul | (Yp + TMoo + a + 5) , 



Q4 =A'Yp5^ ^H,^ +■!<"( Ju + GLn^uZ, ^ +p-m1 ) +panu1, (YpO + (Yp + t:moc)5). 
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Appendix B. Lyapunov functionnal 



In order to obtain a Lyapunov function for system (8|TTi, we approached the problem in a backward 
manner as described in |21, Chap. 4 - p. 120]. We investigated an expression for the derivative <i> 
and went back to chose the parameters of «1> so as to make <!> negative definite. This is useful idea to 
searching for a lyapunov function. Using this method we can derive the global stability in proposition 
[3] The Lyapunov function for the system (8|TTi is: 



4> 



1 (^IP 



el 



1+2 



1 (^ip 



P 



0192 + 9193 + 6293 



1 



9i94 + 29294 



2Tp 
8 



9394, 



where 9i 



,Y« + P^oo+A' ^, 
92 = M — Moo, 93 = p — poc, 94 = — boo, with s\ = max(ri , Tj) such that 



p25«i(l+2i±^) (Yp+A')'(2|;)' [(5+A<)(^;;J^ + l) + (o + 5 + A<)e+2p«„ 



72 = 



1+2^ 
2yp 



'pj V^YP 



e(2 + e] 



1+2 

_5_ 
2Yp 



S+Y„ 8 A 




a 












2 


/ PPoo 






Vy"+p^~ 


+ IJJ 






This Lyapunov function <!> is positive when ^1 +2^^j ^ sf" "a^- ^'■^ derivative along the solutions 
to the system (8pTi is 



P2^7^. 



Y„ + pA„+^ 



9t-pMoo 1+2 



,Y„ + p(A„o + 9i 



9i92 



'2(Y„ + p(Aoe + 9i)+Tp)(Y, + p(Ao„ + 9i)+8) 



~V2Yp 



+ 

Y» + pAoo+A' 

Yp 2Yp 



1 ) +(o + 8 + A()^+2pM. 



Yp 03-8 



e4-(Yp+A') 



_5_ 
2Yp 



Y„ + p(A„„ + 9i) 

Wp 



9i 



6164 



9i93. 



<i> is non-positive. Furthermore, <i> = if and only if 9i = 92 = 93 = 94 = 0. LaSalle's invariant principle 
then implies that the unique equilibrium is globally asymptotically stable in the stable subset defined in 
([I3JII22). 
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